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Actomyosin actively generates contractile forces that provide the plasma membrane with the 
deformation stresses essential to carry out biological processes. Although the contractile property 
of purified actomyosin has been extensively studied, to understand the physical contribution of 
the actiomyosin contractile force on a deformable membrane is still a challenging problem and of 
great interest in the field of biophysics. Here, we reconstituted a model system with a cell-sized 
deformable interface that exhibits anomalous curvature dependent wrinkling caused by actomyosin 
cortex underneath the spherical closed interface. Through the shape analysis of the wrinkling 
deformation, we found that the dominant contributor on the wrinkled shape changes from bending 
elasticity to stretching elasticity of the reconstituted cortex by increasing the droplet curvature 
radius of the order of the cell-size, i.e., tens of micrometer. The observed curvature dependence was 
explained by the theoretical description of the cortex elasticity and contractility. Our present results 
provide a fundamental insight on the deformation of a curved membrane induced by the actomyosin 
cortex. 


I. INTRODUCTION 

Many types of eukaryotic cells have a thin shell-like 
structure, called the cell cortex, that underlies the cell 
membrane and consists of actin filaments (F-actin) and 
type II myosin. Whereas the static stiffness of the ac¬ 
tomyosin cortex maintains the cellular shape against ex¬ 
ternal stresses, the active contractile forces of the cor¬ 
tex play crucial roles in a variety of biological processes, 
such as cell motility [T], cell division [2], embryonic de¬ 
velopment [3], wound healing [4] and cancer metasta¬ 
sis [5]. The actomyosin cortex is interconnected with 
the plasma membrane by actin-related proteins, such 
as ezrin-radixin-moesin (ERM) proteins [6]. However, 
the mechanism of the consequent membrane deformation 
caused by force transduction between the cell cortex and 
cell membrane remains unclear because the membrane- 
cortex complex contains highly intricate components and 
connections [ 7 ]- To understand the intricate roles of the 
actomyosin cortex in cellular processes, bottom-up or re¬ 
constitution approaches with a small number of essen¬ 
tial components have been adopted. For example, actin 
polymerization and depolymerization with actin-related 
proteins in purified systems [8] , the rheological measure¬ 
ments of purified F-actin solution [9], the motor activity 
of myosin from the perspective of its interaction with 
actin ini, and single molecule analysis m have been 
examined thus far. As reconstitutions of biomimetic ar¬ 
tificial systems, the surfaces of giant unilamellar vesicles 
[IMS] or beads [m ng have been utilized frequently 
to investigate the physicochemical properties of interfa¬ 
cial deformation or symmetry breaking of the surround¬ 
ing actin cortices. On the other hand, actomyosin en¬ 
capsulated within a cell-sized confinement exhibits dis- 
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tinct structures, such as network, cortex, and aster for¬ 
mations [13 nnig. However, the fundamental issues 
regarding the actomyosin cortex, i.e., how the mechan¬ 
ical contractile forces exerted by the actomyosin affect 
the deformable plasma membrane and its shape, remain 
poorly understood because of the difficulties of in vitro 
reconstruction, such as those posed by the establishment 
of effective coupling between actomyosin and membranes 
[2Qll22] . Therefore, we sought to develop a model sys¬ 
tem that can exhibit the primitive aspects of cell mem¬ 
brane deformation due to the contractile force generated 
by membrane-coupled actomyosin. 

II. MATERIALS AND METHODS 
A. Cell culture 

Mass culture of Amoeba proteus was performed as de¬ 
scribed previously [23]. In brief. Amoeba proteus was cul¬ 
tured in KCM medium (7 mg/1 KCl, 8 mg/1 CaCl 2 , and 
8 mg/1 MgS 04 ' 7 H 20 ) at 25°C and fed with Tetrahymena 
pyriformis. To avoid contamination with T. pyriformis 
cells were starved for at least 3 days before use. 


B. Preparation of the actomyosin fraction 

The actomyosin fraction was prepared by modifying 
a previously described method [24]. All procedures 
were performed at 2°C. A. proteus (10 g) was sus¬ 
pended in a cell washing solution (2 mM 0,0'-Bis(2- 
aminoethyl)ethyleneglycol-N,N,N' ,N'-tetraacetic acid 
(EGTA), 2 mM MgCl 2 , and 20 mM Piperazine-1,4- 
bis(2-ethanesulphonic acid) (PIPES)-KOH, pH 7.0) 
and centrifuged at 6,000for 2 min, and then the 
supernatant was removed. The resultant precipitate was 
centrifuged at 600,000for 20 min to obtain an actin 
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rich solution. The precipitate was suspended in a 3 M 
KCl solution (3 M KCl, 2 mM MgCl 2 , 1 mM threo-1,4,- 
dimercapto-2,3-butanediol (DTT), 20 /ig/ml leupeptin, 
20 /ig/ml pepstatin A, 0.2 mM phenylmethylsulphonyl 
fluoride (PMSF), 20 mM imidazole-HCl, pH 7.0) and 
centrifuged at 400,000for 10 min. The resultant 
supernatant was dialyzed against 50 mM KCl, 2 mM 
EGTA, 2 mM MgCl 2 , 1 mM DTT, 0.2 mM PMSF, 20 
mM imidazole-HCl, pH 7.0 for 5 h. The dialyzed solution 
was centrifuged at 20,000 for 5 min, and the super¬ 
natant was discarded. The precipitate was suspended in 
a 300 /il EMPA solution (5 mM EGTA, 6 mM MgCl 2 , 
1 mM DTT, 2 mM adenosine-5'-triphosphate (ATP), 
30 mM PIPES-KOH, pH 7.0), and 30 /il actin rich 
solution was added. The obtained fraction dissolving 
the actomyosin contains 4.5 mM EGTA, 5.4 mM MgCl 2 , 
0.9 mM DTT, 1.8 mM ATP and 27 mM PIPES. The 
protein concentration in the actomyosin fraction was 
determined using the method of Bradford |25] with 
bovine serum albumin as the standard, and the ratios 
of actin or myosin to the total protein were measured 
using SDS-PAGE [26] and the Fiji software package 
(http://fiji.se/wiki/index.php/Fiji). Each concentration 
of actin and myosin was calculated relative to the 
concentration of total protein and the ratio of actin 
or myosin to the total protein. Immunoblotting was 
performed according to the standard protocol with 
anti-/3-Actin (Poly6221; BioLegend) and anti-Myosin H 
from Amoeba proteus. Horseradish peroxidase (HRP)- 
conjugated secondary antibodies (Poly4053; BioLegend 
and Poly4064; BioLegend) were used for detection. 


C. Prearation and observation of the actomyosin 
droplets 

Gell-sized water-in-oil (W/0) droplets with a lipid 
monolayer were prepared as described previously 
[27[ l28] . In brief, the amphiphile l,2-dioleoyl-3- 
trimethylammonium-propane chloride (DOTAP, Avanti 
Polar Lipids) was dissolved in chloroform at a concen¬ 
tration of 10 mM. The DOTAP solution was dried in 
a nitrogen stream and allowed to settle under vacuum 
overnight. Mineral oil (Nacarai Tesque) was added to 
the dried films to obtain 1 mM DOTAP, which was soni¬ 
cated at 60°G for 60 min, resulting in dispersed DOTAP 
in oil. Finally, to obtain cell-sized W/0 droplets, 3% 
(vol./vol.) actomyosin fraction was added to the oil solu¬ 
tion, and emulsification was performed by vortexing and 
pipetting. To obtain transmitted light images the pre¬ 
pared droplets were observed using a microscope (Eclipse 
Ti, Nikon) at room temperature and recorded using an 
sCMOS camera (ORCA-Flash4.0, Hamamatsu). To re¬ 
veal the distribution of actomyosin fraction in the W/0 
droplets, the actin filaments were stained with 10 nM 
Acti-stain 488 phalloidin (Cytoskeleton). Note that the 
phalloidin prevents actin filaments from depolymeriza¬ 
tion. Fluorescent time-lapse images were collected using 


a microscope (1X71, Olympus) equipped with an EM- 
GCD camera (iXon, Andor) and a confocal scanner unit 
(CSU-Xl, Yokogawa). In the experiments using droplets 
in the absence of myosin, a mixture of 300 /il EMPA so¬ 
lution and 30 /il actin rich solution was used in place of 
the actomyosin fraction. Note that a myosin H inhibitor, 
blebbistatin, is ineffective in inhibiting myosin in Amoe- 
bozoa m- The obtained images were processed via Fiji 
and analyzed via custom routine in Igor Pro (WaveMet- 
rics, Oregon). 


III. RESULTS AND DISCUSSION 
A. Experimental results 

In the present study, we fabricated a spherical (in its 
initial state) cell-sized deformable interface composed of 
a lipid monolayer imiiiisni that encapsulates the ac¬ 
tomyosin fraction (Fig. [^a)). As shown in Fig. [^b), 
the actomyosin fraction used here consists of actin and 
myosin H extracted from Amoeba proteus [24|. The defor¬ 
mation stress on the actin network cross-linked by rigor 
state myosins [20| or myosin bipolar thick filaments [24] 
was exerted by the motive force of active myosin motors 
through the consumption of ATP, as it was confirmed by 
the observation that a droplet in the absence of myosin 
molecules did not exhibit any deformation. To establish 
nonspecific interconnection between the actomyosin and 
the interface, we employed electrostatic attraction m 
using a positively charged lipid, DOTAP. In this solu¬ 
tion, the positively charged lipid monolayer attracts neg¬ 
atively charged actin filaments, accompanied with myosin 
thick filaments as an actomyosin complex. In this situa¬ 
tion, we confirmed the suitable experimental conditions 
that actin and myosin form an actomyosin cortex struc¬ 
ture underneath the lipid interface and cause interfacial 
deformation by the actomyosin contraction, as 3 mg/ml 
actin, 6 mg/ml myosin, and pure DOTAP lipid mono- 
layer, by tuning the actin concentration, myosin concen¬ 
tration, and DOTAP fraction [32|. If the myosin con¬ 
centration is higher than 6 mg/ml, e.g., 8 mg/ml, the 
actomyosin detached from the lipid layer and condensed 
to form an aggregate that was independent of the droplet 
boundary, as was reported in previous studies [iniET]. 

The interface of the actomyosin droplet was clearly de¬ 
formed under the above-confirmed appropriate condition 
(Fig-B c)), which should be noted that the energy cost 
to deform the interface of a micrometer sized water-in-oil 
droplet is lO^-IO^ times higher than for an elastic lipid 
membrane [33H35] . Nevertheless, it is obvious that the 
droplet oil/water interface deforms because the wrinkles 
revealed by the contrast in the transmitted light defi¬ 
nitely corresponds to the deviation of the lipid layer at 
the oil/water interface from the initial position. This fact 
was further proved by the fluorescence microscopy with 
different dyes for the lipid interface and actomyosin as 
shown in Fig. (see below for the details). For a typical 
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FIG. 1. (Color online) (a) Schematic illustration of a cell¬ 
sized lipid monolayer encapsulating the actomyosin fraction, 
(b) Characterization of the actomyosin fraction; error bars 
represent standard error of the mean (n = 5). (c) Time de¬ 
velopment of the lipid interface induced by the contraction of 
the conhned actomyosin fraction. The wrinkle grew on the 
interface. The scale bar corresponds to 40 /xm. 


case shown in Fig. Ic, initially-spherical droplet exhib¬ 
ited interfacial wrinkling. The inter facial wrinkling was 
not accompanied by apparent symmetry breaking, i.e., 
the strong localization of actomyosin distribution, as oth¬ 
erwise observed in both actin gels around beads [niiiH] 
and reconstituted |2T1 [22] or cellular mill actomyosin, 
owing to the strong adhesion between the actomyosin and 
the lipid interface in the present experimental condition. 


Figure [^a) shows a time development of the wrin¬ 
kle formation visualized by fluorescently-labeled lipid in¬ 
terface and actin with confocal fluorescence microscopy. 
The actin distribution exhibits that the actomyosin inside 
the droplet forms cortex structure underlying the lipid in¬ 
terface. The same fluorescence distributions of the lipid 
interface and actomyosin cortex structure as shown in 
the merged images (Fig. [^a), bottom) and its schematic 
illustrations (Fig. [^b)) clearly verify that the interfacial 
wrinkling was directly induced by the underlying acto¬ 
myosin cortex, which yields lateral contractile force. The 
comparison between the original spherical interface at 
the beginning of the wrinkle formation At = 0 s and the 
wrinkled interface at At = 600 s revealed the apparent 
dents of the deformed part and slight ourward expan¬ 
sion of the entire region at the same time, indicating the 
volume conservation due to the incompressibility of the 
inner aqueous solution (Fig. He)). 
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FIG. 2. (Color online) (a) Confocal fluorescence images 
of lipid interface (labeled with ATTO 488-DOPE) and actin 
filaments (labeled with Rhodamine phalloidin) during the in¬ 
terfacial wrinkling. At = 0 is dehned as the time point at 
which the wrinkling began. Actin hlaments spontaneously 
form a cortex structure and contract underneath the lipid 
interface, (b) Schematic illustrations corresponding to the 
confocal images shown in (a). Red and orange colors denote 
the lipid interface and actomyosin fraction, respectively, (c) 
Comparison between the original spherical interface and the 
wrinkled interface revealed by ATTO 488-DOPE. The original 
interface (At = Os, red) is overlaid on the wrinkled interface 
(At = 600s, cyan). The deformation leads to the slight out¬ 
ward expansion of other part due to the volume conservation 
of aqueous phase. All scale bars correspond to 10 /xm. 


B. Characterization of actomyosin-induced 
deformation 


To quantify and characterize the intrinsic propeties of 
actomyosin-induced deformation, we analyzed a wrinkled 
shape as a perturbative shape from the initial spherical 






















4 


shape. Figure [^a) shows a typical example of the confo- 
cal fluorescence image of an actomyosin cortex. Based on 
the droplet equator revealed here, we calculated the ra¬ 
dial position X of the deformed cortex, which was taken 
from the two dimensional droplet center of mass, as a 
function of circumference length L ranging from 0 to 
27 rR. The mean radius R = (x) angle was defined as 
the rotationally-averaged radial length. The edge detec¬ 
tion was performed by Gaussian fittings of 1024 discrete 
radial intensity profiles in the vicinity of the cortex, in 
which the radial positions x{L) were determined as the 
peak position of each determined Gaussian function. The 
red line depicted on Fig. |^a) (right) indicates the de¬ 
tected shape from the original image (left) via this pro¬ 
cedure. Figure [^b) shows the graph of radial length x 
versus circumference length L, i.e., shape profile of the 
deformation. From the shape profile, we characterized 
the actomyosin-induced deformation by calculating spa¬ 
tial power spectrum and autocorrelation. Gray mark¬ 
ers shown in Figure |^c) are the calculated power spec¬ 
trum x{nAL) — {x{nAL))e‘^'^^^^'^^^ , where 

N = 1024, and nAL = 27 rRn/N, are the sampling 
number, spatial frequency, and circumference length, re¬ 
spectively. The black solid line is the power law fit¬ 
ting, in which the exponent shown as the slope in log- 
log plot reflects the intrinsic properties of the deformed 
shape [33l |36]. Figure [^d) shows the autocorrelation 
{x{Lq)x{Lq -f nAL)), where () denotes the average over 
arbitrary initial position Lq and displacement nAL. Be¬ 
cause of the periodic boundary condition of the shape 
profile (Fig. [^b)), the autocorrelation becomes negative 
at a certain displacement L = nAL. We can define the 
correlation length ^ as the smallest length of L that sat¬ 
isfies autocorrelation = 0, so that ^ represents the char¬ 
acteristic width of a concave in the wrinkled shape m 
For example, in the case of sine function as a deviation 
profile from its mean value, the half-length of a convex 
part, i.e., A/4 (A; the wavelength of the sine function), 
appears as ^ in the autocorrelation. This characteristic 
length of the wrinkles also reflects the underlying mech¬ 
anism of the actomyosin-induced deformation. 


C. Curvature dependence of wrinkled shape 

We examined the spatial properties of the wrinkle 
emergence, which we hypothesize was caused by the lat¬ 
eral contractile stress generated by the actomyosin cortex 
underlying the lipid interface. Note that the distribution 
of actin observed by fluorescence was laterally homoge¬ 
neous throughout the wrinkle development, except for 
the effect of focal depth. Here, we focused on the cur¬ 
vature dependence of the wrinkling deformation induced 
by actomyosin cortex. Along the above described pro¬ 
cedures, the spatial power spectra were calculated for 
various sizes of droplets. Figure [^a) shows the spatial 
power spectra of three typical droplets with various radii. 


Their power-law scaling exponents of spatial frequency q 
change from —4 for smaller droplets (i? ^ 5 /im) to —2 
for larger droplets (R > 20/im). Figure Bb) exhibits 
the convergence of the scaling exponents at —2 with the 
increase in droplet size, ranging from R = 1.9/im to 
53.4 jam (Fig. |^b)). 

These results indicate that the dominant contributor to 
the deformation exhibits transitional behavior in a size 
dependent manner. The crossover size for the appear¬ 
ance of this curvature effect on actomyosin contraction 
is on the order of 10 jam. The exponent —2 of a spa¬ 
tial power spectrum is generally observed in the case of 
interfacial tension-dominated deformation of a planar in¬ 
terface, whereas the exponents —4 and —3 are observed 
in the bending energy-dominated deformation of a planar 
or spherical interface [33l[36]. Therefore, the exponents 
indicate that the deformed shapes of smaller and larger 
droplets are predominantly determined by the bending 
elasticity and contractile interfacial tension of the acto¬ 
myosin shell, respectively. Here, the width of the concave 
interfaces should reflect the elastic and contractile prop¬ 
erties of the actomyosin interacting with the cell-sized 
deformable interface. To gain further understanding of 
the non-specific interaction between the contractile ac¬ 
tomyosin cortex and the deformable interface, we devel¬ 
oped a theoretical model for the deformation based on 
the elastic energy of the actomyosin cortex. Comparison 
between the experimental and theoretical results should 
provide insight into the elastic, contractile, and geometric 
properties of the microscopic actomyosin cortex. 


D. Theoretical model for the onset of wrinkling 

Here, we assume that actomyosin forms a homoge¬ 
neous thin elastic shell-like structure with thickness h 
and Young’s modulus E underneath the surface of the 
droplet with radius R. We also assume that a part 
of the spherical interface is slightly distorted from its 
original spherical shape, with the depth of distortion 
given as 5 (Fig. [^a)). The dent region stores two 
types of elastic energy per unit area: bending energy 
/bend ^ Eh^5 / dx‘^)‘^ ^ Eh^5‘^/v^ and stretching en¬ 
ergy /stretch ^ Ehul^ ^ Ehb^ /R?, where Uab is a com¬ 
ponent of the strain tensor of the shell. In the following, 
we consider the onset of the deformation from the initial 
spherical shape. The elastic energies over the dent re¬ 
gion are obtained as Fbend ^ ^^/bend ^ Eh^6‘^/r‘^ and 
-^stretch ~ J’^/stretch ~ EhS'^r^/E?. Additionally, we 
consider the energy cost when the other part stretches 
outward simultaneously with the partial dent due to 
the incompressibility of the inner aqueous solution (see 
Fig-B c)). When the local interfacial indentation oc¬ 
curs as shown in Fig. [^a), the volume corresponding 
to the dented part is compensated by the slight out¬ 
ward expansion of the other part of the interfacial shell. 
Approximating that the stretching of the undeformed 
part compensates for the dent volume AH ^ 8r‘^, the 
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FIG. 3. (Color online) Detailed analytical procedure of the spatial power spectrum and autocorrelation for the deformed 
equator shown by fluorescently-labeled actin distribution, (a) Typical example of the original fluorescence image (left) and 
the detected edge (right, red line), radial length x is defined as the distance between two dimensional center of mass and the 
detected edge as a function of circumference length L ranging from 0 to 27ri?(/im). Scale bar is 10/xm. (b) Radial length x 
plotted versus circumference length L. (c) Spatial power spectrum and (d) autocorrelation calculated from the shape profile 
shown in ((b)). Correlation length ^ is defined as the smallest length that autocorrelation = 0. 


normal stretching length is Ai? ^ dr‘^ j^ resulting in 
the additional energy cost for the distortion: Fcomp ^ 
Eh{AR/R)^ • i?2(l + Vl - (r/i?)2) - EU^r^/R^ by ne¬ 
glecting the term of 0{{rlR)^). This energy, up to the 
order of (r/i?)^, is smaller compared to Fgtretch by a fac¬ 
tor of , but should not be neglected because the ex¬ 

perimental ratio ^/i?, which corresponds to r/i?, shown 
in Fig. [^b) is not so much small, ranging from 0.15 to 
0.83. 

The interfacial deformation requires a driving force 
along the normal direction to the interface. In our exper¬ 
imental system, the driving force can be the tangential 
contraction provided by the contractile actomyosin shell 
underlying a curved interface Eum]. There is a gen¬ 
eral relationship between the normal stress arr and the 
tangential stress on a curved surface: 

where r' is the radial coordinate [371439] . By integrat¬ 
ing over the shell thickness r' = R ^ i? + /i, the 
normal stress of the actomyosin shell (per unit area) 
can be described with the tangential stress as arr = 
2lR‘^ dr'r'a±j_ ^ 2ha_\__\_/R, where the same ex¬ 


pression is used for the Laplace pressure. Here, the 
tangential stress is assumed to be a constant over 
the shell thickness h. This pressure, which originates 
from the tangential contraction, causes the interfacial 
distortion. The work done by the contraction leads to 
bkcontraction ^ (J7-7-AH ^ hajR. 

Thus, we obtain the total energy cost for the deforma¬ 
tion: 

-^total — -^bend T -^stretch T F comp bFcontraction 

Eh^6‘^ EhS‘^r‘^ / \ ha±±Sr‘^ 

^ i?2 V R ' 

Let us consider the onset of the deformation. In our ex¬ 
periment, the observed dent radii maintained the same 
order of magnitude from the onset time to the end of 
the wrinkling behavior (Fig. |^a)). Therefore, the initial 
dent radii can be regarded as the experimentally observed 
dent radii in the wrinkling period. Immediately after the 
droplet is prepared (t = 0 min), the tangential stress is 
= 0) = 0. As time passes, the stress grows in¬ 
side the shell; consequently, although the energy always 
permits a small agitation of the interface on the order 
of /i, the observed macroscopic indentation of the inter¬ 
face begins at a certain time, t = tc > 0. At this time. 
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FIG. 4. (Color online) (a) Spatial power spectra obtained 
from the typical equatorial shapes of various sized droplets 
represented by pale red, green and blue markers. Correspond¬ 
ing red, green and blue lines are the power hts of the exper¬ 
imental spectra, (b) Values of the power-law scaling expo¬ 
nents in the power spectra on spatial frequency with respect 
to droplet size. The exponents vary from —4 to —2, in corre¬ 
lation with the droplet size. 
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FIG. 5. (Color online) (a) Schematic illustration of the 
situation considered in the model. A portion of the spherical 
interface is slightly dented, where S is the depth of the dent, 
r is the radius of the dent region, R is the initial curvature 
radius of the droplet, and h is the thickness of the actomyosin 
shell underlying the droplet interface, (b) Correlation length 
^ of the wrinkle shape profile along the circumference, plotted 
versus droplet size. Each correlation length is calculated at 
the hnal state of wrinkle development. Red and blue lines 
denote theoretical correlations between the dented radius and 
the droplet radius described in the main text. 


the signature of Ftotai becomes negative in the range of 


0 < r < R. According to these situations, the critical 
conditions for the onset of the dent formation require 
^totai = 0 dFtots,\/dr = 0. These conditions produce 
the following correlation between the dent radius and the 
droplet curvature radius: 

r = {2h^)^Ri. ( 1 ) 

Here, we verified the corresponding experimental values, 
i.e., spatial correlation length ^ of the deformed equato¬ 
rial shapes for various sizes of droplets. ^ was defined 
as the circumferential length at which the autocorrela¬ 
tion first becomes 0, and thus the length ^ denotes the 
width of the concave structure along the circumference. 
Consequently, there is a strong positive correlation be¬ 
tween the correlation length ^ and the droplet radius 
R (Fig. I^b)). The theoretical curvature-radius depen¬ 
dence, indicated by Eq. 0 is consistent with the ex¬ 
perimental slope in the small-radius region with the fit¬ 
ting parameter h = 218 nm, which is comparable to the 
observable results (< 500 nm in smaller droplets); fur¬ 
thermore, this value is close to the cortical thickness of 
living cells j4QH42] and previously reconstituted systems 
j38l [43l l44] . At the same time, the curve Eq. Q de¬ 
viates from the experimental data in the larger-radius 
region of > 10 /im. To elucidate the factors that in¬ 
fluence the deviation, we return our discussion to the 
circumferential power spectra of the deformed shapes for 
various droplet sizes (Fig. [^b)). The scaling exponents 
—4 and —3 indicate that the bending elasticity Fbend 
dominates the deformed shape for < 10 /im, whereas 
the transient distribution and convergence to the expo¬ 
nent —2 indicate that the energies for area stretch (lat¬ 
eral tension), Fgtretch and Fcomp, dominate the defor¬ 
mation for 10 /im < i? < 20 /im and i? > 20 /im. Be¬ 
cause the ratio of the bending term and tension terms 
is idDend/i^stretch = Or Fbend/^d^omp = 

-^bend -^stretch and -^bend -^comp are realized in the 
large-radius region [R > 20/im), where r and R are of 
the same order in experimental values (Fig.j^a)). Thus, 
^bend IS negligible compared to the tension terms for the 
large-radius region. Minimizing the energies in the ab¬ 
sence of the bending energy, we in turn obtain the fol¬ 
lowing equation for the deformation of larger droplets: 


T 



1 

m' 


( 2 ) 


By fitting in the region of i? > 20 /im (blue line in 
Fig. I^b)), we find a±_\_l2E6 ^ 0.023 /im“^. Further¬ 
more, because S is of the order of 1 /im in our experiment 
for R > 20 /im, the above relation results in the quan¬ 
tification of the ratio between the active force generation 
by actomyosin contractility and static elasticity 
i.e., cr^^jE ^ 0.01. Note that Eq. 0 is only verified un¬ 
der the necessary condition a_\__\_l2ES — ll2R > 0, so that 
the energy cost, Ftotai — ^bend? bas a minimum value. In 
the case of a_\__\_l2ES ^ 0.023 /im“^, this necessary con¬ 
dition results in > 20 /im, which is consistent with the 
situation of the larger droplets considered here. 














7 


Although the wrinkles further developed after the on¬ 
set of the deformation and the dent curvature eventu¬ 
ally became negative, the typical pattern of the wrinkled 
shape was determined in the early stage of the inden¬ 
tation. In the late stage of the deformation, the inner 
solution of the droplets evaporated or was ejected (data 
not shown) and the droplets became a bit smaller to re¬ 
lease a part of the pressure increase as in Fig. Be), re- 
suiting in the breaking of the volume conservation inside 
the droplet. In addition to the contractile force of the 
actomyosin shell, such increase in area-to-volume ratio 
in the late stage could induce complex buckling of the 
elastic shell as well. While the wrinkled width ^ or r was 
determined in the early stage considered in the model, to 
investigate the complex deformation mechanism in the 
late stage in detail would be a future work for this sys¬ 
tem. 


IV. CONCLUSIONS AND PROSPECTS 


Reconstituted actomyosin cortex underlying a de¬ 
formable lipid interface successfully caused the wrinkling 
deformation of the interface for the first time. Further¬ 
more, the shape analysis and theoretical descriptions re¬ 
vealed the contributions of the bending elasticity, stretch¬ 
ing elasticity, and contractility of the reconstituted cortex 
structure on the interfacial deformation, which anoma¬ 
lously depends on the cell-sized droplet radius, i.e., the 
interfacial curvature. In other words, we demonstrated 
the physical basis that even the primitive geometrical 
coupling between the deformable interface and confined 
actomyosin exhibits rich variety of interfacial deforma- 
bility. The geometrical origin in our findings should be 
relevant and taken into consideration for not only the 
purified systems but also living systems such as the con¬ 
traction of the cell membrane with the contractile cell 
cortex at the posterior during the amoeboid motion of a 
cell. To study the applicability of the present curvature 
dependence of the contractile cortex in living cells such 
as Amoeba proteus would be one of the most interesting 
future works. With respect to the objective of recon¬ 
stituting cellular motility, the spatiotemporal regulation 
of the structural formations should be addressed as the 
next step of this work by exploring the implementation 
of the frameworks of actin-related molecules and their 
signaling networks. Our present findings on how acto¬ 
myosin behaves within a cell-sized compartment with a 
deformable interface provides an experimental basis for 
this purpose, and it takes an important step toward ar¬ 
tificial cellular deformation in the step-wise reconstitu¬ 
tion of the living state. Together, these investigations 
may lead to the further clarification of complex biologi¬ 
cal phenomena involving geometrically coupled actin and 
myosin functions. 
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Appendix: EFFECT OF INTERFACIAL TENSION 

In the main text, we only considered the elasticity and 
contractility of actomyosin shell. However, the surface 
tension of the lipid interface may also contribute to the 
deformation energy cost, because the lipid interface was 
also deformed to the same shape with the shell, as shown 
in Figure a). In the same theoretical framework used 
in the main text, the contribution of surface tension 7 
for the dent region can be described as the work done by 
Laplace pressure 


‘^Av ~ 

R ~ R 


(A.l) 


where AK ^ 7 rJr^ is the decreased volume by the dent 
defined as negative. On the other hand, the other region 
should slightly expand outward due to the imeompress- 
ibility of the inner aqueous phase, resulting in the radial 
expansion AR ^ ^. This contributes to the interfacial 
energy cost: 


7 AA 


87rj6r‘^ 
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(A. 2 ) 


where AA An {{R^ AR)‘^ — is the approximated 
surface area difference of the droplet before and after 
the outward expansion. Therefore, the interfacial energy 


cost of the deformation becomes Fint — — 27rj^r ^ 

The total energy cost including the elasticity and 
contractility of actomyosin shell is then obtained by re¬ 
placing the lateral stress a_\_± by effective lateral stress 


Eota.l — Fhe.nd -Fstretch H“ F, 




W. 




4- P. 




Eh^5^ 


EhS^r^ 

i ?2 



h5r^ ( 7 \ 

R V h) ■ 


The same fitting in the main text reads ^ ^ 0.01. 
Thus, the experimentally observed wrinkling is realized 
in the case that the contraction exceeds the interfacial 
tension, namely Our experimental result in¬ 

dicates the lower limit of the ratio between the active 
contractility and static elasticity E as > 0.01. 
Note that h ^ 200 nm and 7 ^ lOmN/m result in 
^ ^ 1 nN /m^. 
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